AMENDMENT UNDER 37 C.F.R. § 1.1 1 1 
U.S. Appln. No. 10/766,011 

AMENDMENTS TO THE CLAIMS 

This listing of claims will replace all prior versions and listings of claims in the 
application: 

LISTING OF CLAIMS: 

1 . (currently amended): A gene expression state estimating system for estimating the 

probability of gene expression in each channel, the system including an input device for sending 
gene expression level data, a program-controlled data analyzer, and an output device, wherein 
said data analyzer comprises 

distributed parameter estimating means for estimating distributed parameters of a mixed 
normal distribution shown in the following equation (25) using the gene expression level data 
from said input device, and sending the estimated distributed parameters: 



normal distribution with average jj^O and variance <£*al , <j> (w - //, | a] ) re presents the density 
function of a one-dimensional normal distribution with average u . and variance q^. 
(// 0 ,o- 0 2 )and Cu l5 <y, 2 ) (//u,^ 2 ) are average and variance parameters of first and second 
components, respectively, u represents data points corresponding to sum of gene expression level 



(1 - I) ip {u - p 0 | Oo) + 5<p (u - ji x I oj) 



(25) 



where Mu-fi, | a 2 ,) 1 ^^ 



represents the density function of a one-dimensional 




.and 4 is the mixing ratio, with the assumption that 



Fo < Pi, °o > °* *>t > °> 0 < 5 < 1 



is satisfied, 
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mixing ratio parameter estimating means for estimating a mixing ratio parameter of the 
mixed normal distribution using the gene expression level data sent from said input device and 
the distributed parameters sent from said distributed parameter estimating means, and sending 
the estimated mixing ratio parameter, and 

posterior probability calculating means for calculating the posterior probability of the 
expression state of each gene in each channel using the gene expression level data, the estimated 
distributed parameters and mixing ratio parameter, and sending the calculated posterior 
probability to said output device. 

2. (original): The system according to claim 1 wherein said distributed parameter 

(a 2 a 2 ) 

estimating means estimates the mixing ratio (!;), average (u 0 , Uj), and variance 1 °' 11 by 
applying the mixed normal distribution of two components to data on the sum of the amounts of 
expression of genes located in a region where the difference of gene expression levels X and Y 
of two channels is near 0. 

3. (currently amended): The system according to claim 2 wherein when the median 
value of the absolute difference ]vj|(i = l,...,n) of the gene expression levels X and Y is c M , the 
data on the amounts of gene expression is shown by {uj | j Vj | < c M , .i~l..,..nl . and v represents 
data points corresponding to difference of the gene expression level of the two channels . 
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4. (original): The system according to claim 3 wherein said distributed parameter 

estimating means performs estimation by the use of the estimated ^ * 0# ° x according 

to the following equations (26), (27), (28), and (29) to estimate p ' °*' a P' * ; 
F = <£i - Po ) / 2 (26) 

ii y n i<eN 0 ^7) 

«2 _ 1 «2 1 -2 

4 2 (28) 

where No denotes an index set of data values that satisfies * € ^ I u± < ^ and 
II N °J1 denotes the number of elements. 

5, (currently amended): The system according to claim 4 wherein said mixing ratio 
parameter estimating means estimates the mixing ratio parameter p=(poo, pio, Poi, pn) (where p 0 o 
denotes a mixing ratio when no gene is being expressed in both cells 1 and 2, pi j denotes a 
mixing ratio when any gene is being expressed in both cells 1 and 2, p w denotes a mixing ratio 
when a gene is being expressed in cell 1 but not in cell 2, and q^Sm- represents a mixing ratio 

when a gene is being expressed in cell 2 but not in cell 1) using 9 _ * given 

from said distributed parameter estimating means by applying a two-variable mixed normal 
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distribution of four components shown in the following equation (30) to the-gene expression 
level data {(uj, Vj)|i-1, n} sent from said input device wherein v represents data points 
corresponding to difference of the gene expression level of the two channels, 

Poo3oof u ' v I ®> t PioSio< u ' v | §) + Poigoa.tW' v | 9} + p u g n {u, v | §) 

* p 00 <p(u I 4S| + 2o|)tp(v | 2al) + p 10 <P 2 <» ~ V - p | S 10 ) 

+ Poi<P 2 tu - P> V + V I Sol + Pll«p(» - 2£ I 4p 2 (e x2 - 1) 

+ 46= + 2d 2 e mv | 2dl) (3Q) 

where the above equation satisfies the relationships shown in the following equations 
(31) and (32). 



jT(e A - 1) 



- 1) 
1) + ) 



(31) 



£ 2 (e X - 1) + 4o*. + 2of 



1) 



(32) 



6. (currently amended): The system according to claim 5 wherein said posterior 
probability calculating means calculates the posterior probability of expression of any gene in 
cell 1 and cell 2 for each pair (u, v) of the gene expression level data sent from said input device 

according to the following equation (33) (where f tu * v ! ® } is given by the following 
equation (34)) and the following equation (35) (where gene expression level data are shown in 
the equations (X = x L a + (3 + eu Y - + 6 + sA andf whefe n, % 2 take either 1 or 0, which 
represents the presence or absence of true gene expression in each cell). 
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J?r (t * 1 I p, 9) = PiQgio^ v 1 ft + Pirfii^ v I ®> 

1 ' f to v | p, 9) (3 3) 

f {u, v | p, §) = £ p jk g 3k (u, v | 9) 

(j,k)e{0,l} 2 (34) 



Pr U 2 = 1 | p, 9) = PoigPl^ v 1 §) + Pligil^ y I e> 
1 f (u# v | p, el 



7. (currently amended): The system according to claim 5 wherein said posterior 
probability calculating means calculates the posterior probability indicating an event of 
differential expression between cell 1 and cell 2 according to the following equation (36) 

f (u, v I p f 9) (36 j 
where gene expression level data are shown in the equations (X - T i a + |3 + Su Y - t?q + 
p + eA and ti. T2 take either 1 or 0, which represents the presence or absence of true gene 
expression in each cell. 

8. (original): The system according to claim 6 or 7 wherein said posterior probability 
calculating means judges whether calculations of posterior probabilities of gene expression have 
been made for all the pairs (u, v) of the gene expression level data, and when all the calculations 
have been completed, it ends the process, while when all the calculations have not been 
completed yet, it calculates the posterior probability related to the next gene, such that the 
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calculated posterior probabilities of gene expression in each channel are sent to said output 
device, and 

said output device displays the posterior probabilities of gene expression in each channel. 

9. (original): A gene expression state estimating method for estimating the probability of 
gene expression in each channel based on gene expression level data, comprising the steps of: 

estimating distributed parameters of a mixed normal distribution shown in the following 
equation (37) using the gene expression level data and sending the estimated distributed 
parameters: 

(1 - I) cp (u - no J oj) + |<P (u - Pi j oj) (37) 
where <P ( * I ° ) represents the density function of a one-dimensional normal 

distribution with average 0 and variance a 2 , * Po ' 9 ° 1 and are average mid variance 

parameters of first and second components, respectively, and % is the mixing' ratio, with the 

assumption that *° < ^ °° > 0, 4 > 0, 0 < \< 1 is satisfied> 

estimating a mixing ratio parameter of the mixed normal distribution using the gene 
expression level data and the estimated distributed parameters, and sending the estimated mixing 
ratio parameter, and 

calculating the posterior probability of the expression state of each gene in each channel 
using the gene expression level data, the estimated distributed parameters, and the estimated 
mixing ratio parameter, and sending the calculated posterior probability. 
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10. (original) : The method according to claim 9 wherein said step of estimating the 
distributed parameters further comprises a step of estimating the mixing ratio (§, average (no, 

{ 2 2, 

ji]), and variance °' 1 by applying the mixed normal distribution of two components to data 
on the sum of the amounts of expression of genes located in a region where the difference of 
gene expression levels X and Y of two channels is near 0. 

11. (original): The method according to claim 10 wherein when the median value of the 
absolute difference |vi| (I = l,...,n) of the gene expression levels X and Y is c M , the data on the 
sum of the amounts of expression of genes is shown by {iij | j Vi | < cm, i=l,...,n}. 

12. (original): The method according to claim 1 1 wherein the estimation is performed in 

S a ~2 ~2 

said step of estimating the distributed parameters by the use of the estimated ° 0 ' °i 

according to the following equations (38), (39), (40), and (41) to estimate p ' ° e ' X ; 
$ = % " #o>/ 2 (38) 
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where No denotes an index set of data values that satisfies 1 e ^ I Ui < and 

j|N II 

" "'I denotes the number of elements. 

13. (original): The method according to claim 12 wherein the estimation is performed in 
said step of estimating the mixing ratio parameter p=(poo, pio, Poi, pn) (where p 0 o denotes a 
mixing ratio when no gene is being expressed in both cells 1 and 2, pi i denotes a mixing ratio 
when any gene is being expressed in both ceils 1 and 2, p ]0 denotes a mixing ratio when a gene is 
being expressed in cell 1 but not in cell 2, and goi represents a mixing ratio when a gene is being 

expressed in cell 2 but not in cell 1 ) using ^ ~ °* ' 0 P ' sent from said step of 

estimating the distributed parameters by applying a two-variable mixed normal distribution of 
four components shown in the following equation (42) to the sent gene expression level data {{iij, 
v,))i- 1 n} 

PO-OSoot* v ! §) + Pio9io<^ v I S) + Poigoi(* v | §) + p n g n (u, v I 9) 

- Poo<P(« I *of + 26*)<p(v ) 25g) + Pioq> 2 tu - 0, v - £ | E 10 ) 

+ Poi» 3 -(-U ~ ft v + ja | Eo X + PnqKu - 2p | 4ft 2 (e^ - 1) 

+ 4c 2 , + 2o\) <p (v | 2o 2 ) ^ 

where the above equation satisfies the relationships shown in the following equations 
(43) and (44), 



p 2 (e X - 1} + 40 2 + 2d 2 p 2 (e X - 1) 
[ i* 2 ^ - 1) 



(43) 
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K - 1) + 43§ + 23* 

- - 1) 



(44) 



14. (original): The method according to claim 13 wherein the calculation of the posterior 
probability of expression of any gene in cell 1 and cell 2 is made in said step of calculating the 
posterior probability for each pair (u, v) of the sent gene expression level data according to the 

following equation (45) (where f ^ v ! P' ®* i s given by the following equation (46)) and the 
following equation (47) (where ti, t 2 take either 1 or 0, which represents the presence or absence 
of true gene expression in each cell). 

Pr (t,. - 1 I p, 9} . friQgio^ v 1 §) + fagMfo v J h 

Un, v j p, 9) (45) 

f Cu, v I p, §) = 2 PjkST jk(^ v 1 9) 

(j,k) e {0,i} 2 (46) 

Pr h = 1 I p 9} ■= PQigQi^ v I ® + Pugn^ v I §) 

3 v | p,§) (4?) 

15. (original): The method according to claim 13 wherein the calculation of the posterior 
probability indicating an event of differential expression between cell 1 and cell 2 is made in said 
step of calculating the posterior probability according to the following equation (48) 

Pr^ * t 3 I p, 9) = gMgig^ * 1 §) * Poigoifa* v. i 5) 

f v | p, 9) (48) 
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where i\, t 2 take either 1 or 0, which represents the presence or absence of true gene 

expression in each cell. 

16. (original): The method according to claim 14 or 1 5 wherein it is judged in said step 
of calculating the posterior probability whether c alculations of posterior probabilities of gene 
expression have been made for all the pairs (u, v) of the gene expression level data, and when all 
the calculations have been completed, the process is ended, while when all the calculations have 
not been completed yet, the posterior probability related to the next gene is calculated. 

17. (currently amended): A gene expression state estimating program for estimating the 
probability of gene expression in each channel based on gene expression level data, the program 
instructing a computer to execute the steps of: 

estimating distributed parameters of a mixed normal distribution shown in the following 
equation (49) using the gene expression level data and sending the estimated distributed 
parameters: 

a - i) q> (a - ]i 0 i og) + |cp (u - pi | ol) (49) 

where <f>(u-fj 0 I <?l ) ^ I represents the density function of a one-dimensional 
normal distribution with average u^O and variance <r 0 2 <f> (u - //, | of ) represents the density 

function of a one-dimensional normal distribution with average u^ and variance cr^, ^ 0f 
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and ^i' are average and variance parameters of first and second components, respectively, 
u re presents d ata points corresponding to sum of gene expression l evel from two channels and t, 

is the mixing ratio, with the assumption that Po < ^ °° > °' ai > °' 0 < 5 < 1 is 
satisfied, 

estimating a mixing ratio parameter of the mixed normal distribution using the gene 
expression level data and the estimated distributed parameters, and sending the estimated mixing 
ratio parameter, and 

calculating the posterior probability of the expression state of each gene in each channel 
using the gene expression level data, the estimated distributed parameters, and the estimated 
mixing ratio parameter, and sending the calculated posterior probability. 

18. (currently amended): A computer-readable recording medium storing a gene 
expression state estimating program for estimating the probability of gene expression in each 
channel based on gene expression level data, the program instructing a computer to execute the 
steps of: 

estimating distributed parameters of a mixed normal distribution shown in the following 
equation (50) using the gene expression level data and sending the estimated distributed 
parameters: 

(1 - |) <p (u - p 0 I oj) + (u - Pl I oj) 
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where ^ (u - ju 0 \ <rj) ^* I represents the density function of a one-dimensional 

normal distribution with average u_gQ and variance © J cr 2 s .<p(u- }i x | cr 2 ) -represents the density 

function of a one-dimensional normal distribution with average m and variance a f > 

(Po* °o ) ^ j CJ . } ^ are average and variance parameters of first and second components, 
respectively, u represents data points corresponding to sum of gene expression level from two 
channels and £; is the mixing ratio, with the assumption that 
p 0 < p^o? >0,oj > 0, 0 < % < l issatisfiedj 

estimating a mixing ratio parameter of the mixed normal distribution using the gene 
expression level data and the estimated distributed parameters, and sending the estimated mixing 
ratio parameter, 

calculating the posterior probability of the expression state of each gene in each channel 
using the gene expression level data, the estimated distributed parameters, and the estimated 
mixing ratio parameter, and sending the calculated posterior probability, and 

outputting the calculated posterior probability. 
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